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Abstract. We study theoretically the effects of disorder on Bose-Einstein 
condensates (BEC) of bosonic triplon quasiparticles in doped dimerized quantum 
magnets. The condensation occurs in a strong enough magnetic field, where the 
concentration of bosons in the random potential is sufficient to form the condensate. 
The effect of doping is partly modeled by a 5 - correlated distribution of impurities, 
which (i) leads to a uniform renormalization of the system parameters and (ii) produces 
disorder in the system with renormalized parameters. This approach can explain 
qualitatively the available magnetization data on the Tli-zK^CuCbj compound taken 
as an example. In addition to the magnetization, we found that the speed of the 
Bogoliubov mode has a maximum as a function of x. No evidence of the pure Bosc 
glass phase has been found in the BEC regime. 



PACS numbers: 75.45. +j, 03.75.Hh, 03.75.Kk 



The effects of disorder in dimerized quantum magnets in mean field approximations 3 
1. Introduction 

The effects of disorder on the properties of Bose-Einstein condensates present interesting 
problems both for theoretical and experimental physics [H |2l [31 HI [6]. Disorder 
is important in various systems of real particles such as superfluid 4 He, cold atoms in 
optical lattices, and quasiparticles such as polaritons [7] and excitons [8] • These systems 
are well-suited for experimental studies, however the theory of disordered ensembles of 
interacting bosons is complex and there are essentially no exact solutions even in one 
dimension [9]. To approach this problem, Yukalov and Graham (YG) developed a 
self-consisted stochastic mean field approximation (MFA) [10] for Bose systems with 
arbitrarily strong interparticle repulsion and arbitrary strength of disorder potential. It 
was shown that, in general, the Bose system consists of following coexisting components: 
the condensate fraction, po, the normal fraction p^, the glassy fraction pa, and, in 
addition, can be characterized by the superfluid density p s . In the limit of asymptotically 
weak interactions and disorder the known results, obtained in pioneering work by Huang 
and Meng [11] (HM) are reproduced by the YG theory. An interesting question here 
concerns the problem about the existence of a pure Bose glass (BG) phase, i.e. the 
phase where the condensate fraction is nonzero, while the superfluid fraction is not 
yet present. Note that Ref. [6] introduced an alternative definition of the gapless BG 
phase, having localized short-lived excitations and vanishing superfluid density with 
a continuous transition to the normal phase at finite temperature. Even without 
disorder, the condensate is depleted by particle-particle interactions and temperature. 
The inclusion of random fields depletes the condensate further and, possibly, creates the 
glassy fraction. 

As it was understood recently, a new class of BEC can be provided by spin- 
related quasiparticles in magnetic solids such as intensively pumped magnons [12] or 
triplons in the dimerized quantum magnets in the equilibrium [13] . In the magnets, the 
effect of disorder, which can be produced by admixing other chemical elements, can be 
rather strong to be seen in the physical properties such as the temperature-dependent 
magnetization. The so far most investigated compound showing BEC of triplons is 
TlCuCl 3 . To study the effect of disorder, solid solutions of quantum antiferromagnets 
TICUCI3 and KCuChj, i.e. Tli-^K^CuCls have been experimentally investigated recently 
[HI [151 [16] at low temperatures T. The zero-field ground states of TlCuCl 3 and KCuCl 3 
are spin singlets with excitation gaps A st = 7.1 K and A st = 31.2 K, respectively and 
the magnetic excitations are spin triplets. Triplons arise in magnetic fields H > H c , 
where H c is defined by condition of closing the gap by the Zeeman splitting, that is 
A st = OP'bHc, where g is the electron Lande factor and p,B is the Bohr magneton. In 
the mixture Tli_-K-CuCl3 the induced magnetization M exhibits a cusplike minimum 
at a critical temperature T C (H) for fixed magnetic field H > H c similarly to the parent 
compound, and can be successfully explained in terms of triplon BEC [T7J [18j [19] . 

For a theoretical description it is natural to assume that for weak doping a; < 1 
in the mixed system Tb_a.Ka.Cud3 a small admixture of potassium forms a disorder 
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potential. Consequently, the recently developed theories of "dirty bosons" [T| l_j [TU| [TTj 
can be applied to study the BEC of triplons in Tb^KzCuCls. Here the following 
natural questions arise. For example, what is the correspondence between admixing 
parameter x and the properties of the disorder potential ? What are the experimental 
consequences of the disorder ? Yamada et al. [15] analyzed the electron spin resonance 
spectrum in Tb.ajK^CuCFj and concluded that there is a Bose glass - BEC transition 
near a critical magnetic field. Although this interpretation might need a further analysis 
(see discussion in Refs. [20]) it would be interesting to study the influence of the glassy 
phase, or more exactly, of the glassy fraction pa on the magnetization. Note that even 
the existence of a pure Bose glass phase still is a matter of debate even in theoretical 
approaches. For example, it may be predicted by the approach used by Huang and 
Meng [H] if one extends their formulas from weak disorder to a strong one. On the 
other hand, no pure Bose glass was found in Monte - Carlo simulations [21] for atomic 
gases, but predicted for triplons at T = by Nohadani et al. [22]. 

Here we develop a theory of the disorder effects on the BEC of triplons taking 
Tb_a.Ka.CuO3 as a prototype for studies of specific properties. For example, in atomic 
gases considered in Refs. [T_] [H] , the chemical potential fi is determined self-consistently 
with fixed number of atoms, while in the triplon gas the chemical potential is a given 
external parameter controlled by the applied magnetic field and the number of triplons 
is conserved in the thermodynamic limit. To clarify the terms, we underline that the 
number of magnons may vary but that of triplons may be tuned and kept fixed, which 
makes possible the BEC of the latter. 

The paper is organized as follows. In Sections II and III we outline the YG and 
HM approaches valid only for T < T c and extend it for the triplon system. The shift 
of T c due to disorder and the normal phase properties will be discussed in Section IV. 
Our numerical results will be presented in Section V. Conclusions will summarize the 
results of this work. 



2. Yukalov- Graham approximation for disordered triplons 

In the following we reformulate the Yukalov- Graham approximation to the triplon 
system with arbitrary disorder. The Hamiltonian operator of triplons with contact 
interaction and implemented disorder potential V(r) is given by 

U , , + ^ 



H = d 6 r 



^(r) (k - im + V(r)) ^(r) + - {^(r^r))' 



(1) 



where ip(r) is the bosonic field operator, U is the interparticle interaction strength, and 
K is the kinetic energy operator which defines the bare triplon dispersion e^. Since the 
triplon BEC occurs in solids, we perform integration over the unit cell of the crystal with 
the corresponding momenta defined in the first Brillouin zone. Below we assume that the 
bare spectrum remains coherent in the presence of disorder and consider it as a simple 
isotropic one: = k 2 /2m, where m is the triplon effective mass. The distribution of 
random fields is assumed to be zero - centered, (^(r)) = 0, and the correlation function 
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R(r — r') = (V(r)V(r')). Here and below we adopt the units k B = 1, H = 1, and V = 1 
if not stated otherwise for the unit cell volume. 

To describe Bose condensed system where the global gauge symmetry is broken, 
one employs the Bogoliubov shift: 



V(r) = Vi^W+^iW (2) 

where the condensate density Po( r ) is constant for the homogeneous system, po( r ) = Po- 
Since by the definition the average of ip^(r)ijj(r) is the total number of particles: 



N = d 3 r(^(r)ip(r)) (3) 
Jv 

with the density of triplons per unit cell p = N/V, from the normalization condition 

P = Po + Pi (4) 
one immediately obtains 

\ [ cZ 3 r<v4(r)Vi(r)). (5) 



Therefore the field operator ip\(r) determines the density of uncondensed particles. 

The YG approximation is formulated in representative ensemble formalism, which 
includes two Lagrange multipliers, /i and pi, defined as: 

dpuo dpi 

where Q is the grand thermodynamic potential. It was shown that disorder would not 
change the explicit expressions for chemical potentials, obtained earlier [23] in Hartree- 
Fock-Bogoliubov (HFB) approximation without disorder, 

Po = U{p + Pi Pi = U{p + P\ ~ or), (7) 

1 r « 



where a = — d r (^i(r)tfji(r)) is the anomalous density. The total system chemical 
V Jv 

potential mu, related to the total number of particles as N = —dQ/dfi, is determined 
by 

fip = //ipi + /i po- (8) 

Clearly, when the gauge symmetry is not broken, i.e. po = 0, a = 0, p\ = p, both p 
and pi coincide giving p = pi = 2Up. 

In contrast to homogeneous atomic gases considered in Refs. flDlITT] . where p is fixed 
and p(p) should be calculated as an output parameter, in the trip Ion gas the chemical 
potential is fixed by the external magnetic field, while the density p = p(p) should be 
calculated self consistently. In fact, in a system of triplons p characterizes an additional 
direct contribution to the triplon energy due to the external field H and can be written 
as 

p = gp B H - Am, (9) 
which can be interpreted as a chemical potential of the S z — — 1 triplons. 
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The magnetization is proportional to the triplon density 

M = gn BP (10) 
with p is defined from (jSJ) as 

p = - (/iipi + (jiopo) (11) 
H 

where /i and p,\ are given in © and the densities p , p\ must be calculated self 
consistently. 

It is well known [23] that the disorder field leads to formation of a glassy fraction 
with the density pa- In this approximation each of pi and a are presented as 

Pi = Pn + Pg] (T = o- N + p G (12) 

where p^ and are the normal and anomalous densities without disorder. In the YG 
method, based on HFB approximation, the following explicit relations can be obtained 
HH: 



(Am) 3 / 2 r d s k £k + A 

37T 2 J (27f) £k 



P^v = _ ? + / 7TT3/s(^fcJ— , (13) 



(Am) 3 / 2 [ d 3 k „ . c .l 



with the Bose distribution of Bogoliubov excitations /b(^) = l/(e £fc / T — 1) having the 
dispersion 



^ = vW^fc + 2A. (15) 

For small momentum /c the dispersion is linear, = c/c, and the speed of the Bogoliubov 
mode 

\/A , . 

Vm 

The self energy A is determined formally by the same equation as in the case when 
the disorder is neglected, 

A = U(po + a) = U(p- Pn + o- n ). (17) 

The contribution from the disorder potential is hidden in the density of the glassy 
fraction 

p G = ^ J d 3 r((M*)Mr))) (18) 

where the double angle brackets mean the stochastic average. In general, the calculation 
of pa is rather complicated, however, for the 5 - correlated disorder i.e. for the white 
noise, 

((V(r)V(r'))) = R5(r-r'), (19) 
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equation (ITS!) is simplified as [TU] 

_ Ro(p ~ Pn) , , 

PG ~ R + 7(l-R )^- (20) 

The density of condensed fraction can be found by inserting (12) and (20) into the 
normalization condition fl4]). The result is 

7(1 -Rof/\ P - PN ) 
p0 ~ R + 7(l-R )^ ■ (21) 

In Eqs. (1201) and (1211) we introduced the dimensionless parameter Rq as 

7Rm 2 , . 

Ro = — r =. (22) 
47TV mA 

One can see from Eqs. ( !20|) and ( T2~TT) that the glassy fraction is proportional to the 
condensed one, 

pG ~ 7{i-R yn- (23) 

The system of Eqs. (JTJ) , ([8]) , (JT3l)- (fT9l) are the basic of YG approximation. 

An interesting quantity, crucial for determining the Bose glass phase, is the 
superfluid density, p s . In general it is defined as a partial density appearing as a response 
to a velocity boost 

ft = S57S5^> (24) 

where P v is the total momentum of the system, dependent on the macroscopic velocity 
v. Referring the reader to original papers JTOj [11] we bring below analytical expression 
obtained there for p s in the case of white noise random potential 

Ps =,- — -—, (25) 

i r k 2 d 3 k 

lm~J (27r) 3 sinh 2 (^ fc /2T) 
Note that YG approach is valid for arbitrary strength of the interaction potential 
U, and for arbitrary strong disorder. For the weak interactions it leads to pioneering 
Huang- Meng approach which will be extended to the "dirty triplons" in the next 
section. 



Qn = — I 77^—727^17^- ( 26 ) 



3. Huang-Meng approximation 

For completeness, we present here the results for the Huang-Meng approach, based on 
the so called Hartree Fock Popov (HFP) approximation which has been widely applied 
in the literature to describe the BEC of triplons [171 [19]. The basic equations of this 
approach can be obtained by neglecting the anomalous density a, which leads naturally 
to the single chemical potential p, = fio = P>i- Namely, one finds from (JTJ), (JH]) and ( TT71) 



A = Upo, p = U(p + Pl ). 



(27) 
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From these equations and (1121) one obtains following main equations for the self energy 
A: 

A = p-2U(p N + p G ), (28) 

where p^ is formally given in (I13p . and po is determined by the first equation in fl27|) . 
The glassy fraction can be obtained from (1201) in the linear approximation by R assuming 
weakness of interparticle interaction [10] E] 

where a s = Um/An is the s - wave scattering length. Inserting (|29|) into (|28l) we can 
rewrite the former as 

m 2 R\ r A . . 

A = H-2Up N - n L_ . 30 
27r-v/m 



To evaluate the densities one has to solve nonlinear algebraic equation ( 130]) . where p^ 
is given formally by (!I3|) . with respect to A. Next, by inserting the result into ( |27l) 
and ( 129]) one obtains the density of condensed triplons p an d the glassy fraction p G , 
respectively. The total density can be evaluated then by the normalization condition 
P = Po + Pn + Pg- Equations (1231) . ( 126]) for the superfluid density are formally the same 
in both approximations. 

4. The shift of the critical temperature due to disorder and the T > T c 
regime 

It is well known that the critical temperature of BEC, T c for an ideal gas is given by: 

where p c is the total density of triplons near the critical temperature of BEC for pure 

system, 

p c = n/2U, (32) 

with C{x) being the Riemann function. Equation (132]) directly follows, from Eqs. (J7J), 
(JHJ) or (128]) by setting p N = p and p = p G = 0. 

Clearly, any type of interaction is expected to modify T c . In general, these 
modifications are related to the interparticle interactions as well as to the disorder 
potential. Both approaches, considered here give a zero shift due to the boson-boson 
repulsion. However the shift due to the S— correlated disorder (fT9]) . AT C = T c — T c ° is 
given as [TQl EH] 

AT r 2u 



T c ° 9tt' 
where the dimensionless disorder parameter v 

1 



1/3 T ' 
Pc Moc 



(33) 



(34) 



The effects of disorder in dimerized quantum magnets in mean field approximations 9 



c 



is introduced with the localization length 

Lloc = WR- (35) 
For practical calculations we rewrite T c in Eq. (|33|) . which is in a good agreement with 
perturbative estimates [2] as well as with Monte Carlo simulations [26], as an explicit 
function of effective mass m, the interaction strength U, critical magnetic field H c , 
disorder parameter u, and external field H as follows: 

9tt - 2u ( V2g/2 B (H - H ( " 

Now we proceed to consider the triplon density in the normal state in the T — T c ^> AT, 
temperature range. The dirty bosons in the normal phase where the gauge symmetry 
is not broken, are yet poorly studied. For R = with p = pc = o = the triplon gas 
behaves like an "ideal gas" with an effective chemical potential p e ff , and the density 

d 3 k 1 
(27r) 3 exp((e fc -/i cfr )/T) - 1 
Although p, e ff is not accurately known it depends in general, on p, as well as on R. For 
the pure case MFA [27] gives p e s(R = 0) = p — 2Up. The contribution from the disorder 
potential has been studied neither in YM nor in HM approaches. Therefore, to make 
the calculations self consistently, we have to use 

f d 3 k 1 
P(T >Tc) -J (27f exp ((e k - p - 2U P )/T) - 1 ' (38) 
which yields the density p as a solution of the nonlinear equation (138j) . 



p(T >T C )= I „ Wfm r - (37) 



5. Results and discussions 

In the calculations below, the energies are measured in Kelvin, the mass in K _1 , the 
densities are dimensionless and the Bohr magneton is fiB = 0.671668 K/T. As to the 
strength of disorder potential R, it has units K -2 while the disorder parameter u, defined 
in Eq. (1511) is a number supposed to be less than one, v < 1. As a material parameter, 
we use mean dimer-dimer distance in TlCuChj r^d = 0.79 nm [19]. 

To perform numerical calculations in the YG approximation, assuming that fi, U, 
m, and R are given parameters, we use following strategy, (i) By inserting (J7J, (fl2]) . 
(120]) and (T2T1) into (ITT]) we obtain quadratic algebraic equation with respect to p and 
solve it analytically, (ii) By using this p(fi,R,A) and ( IT3]) . ( IT4l) in ( IT71) we solve the 
latter numerically with respect to A, and (iii) by inserting this A back into p(p, R, A) 
we find the magnetization from ( ITOl and evaluate other densities like po and Pg from 
(|2Q) and d2U>. 

In Figure UJ we present as an example the total triplon density p(T) for a clean and 
strongly disordered [y = 0.45, see Eq.( l34l) ) TlCuCFj, obtained in the YG approximation 
assuming that the total effect of the doping leads only to the randomness in the triplon 
subsystem. 
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I i l i l i l i l i I 

2 4 6 8 10 
r[K] 

Figure 1. The total triplon density as a function of temperature in the YG 
approximation for two values of v. Here the following set of parameters m=0.0204 
K _1 , A st = 7.3 K,U = 313 K, and g = 2.06 [Hj valid for TlCuCl 3 is used. 

The calculation of other quantities using the same assumption shows that the 
disorder leads to a decrease in the condensed and superfluid fractions, thereby increasing 
the glassy one. This tendency is quite natural, since the localization effects prevent 
particles from going into BEC. However, the increase in pc is so weak that along 
with p the total number of triplons p is also decreased with increasing the strength 
of disorder potential R. Bearing in mind that p is proportional to the magnetization M, 
and v is assumed to be approximately proportional to x, and comparing FigJT]with the 
experimental data illustrated in Fig J2] one may conclude that the agreement between the 
theory and the experiment is unsatisfactory since the main features of the experimental 
results are not reproduced there. As it is seen in Figj2]the disorder leads to an increase 
in the magnetization and, hence, in the total triplon density. This is accompanied 
by the decrease in the transition temperature. We therefore conclude that while the 
triplon gas can be considered similarly to atomic gases for which the considered mean- 
field approximations were developed, some further additional specific material - related 
properties of the dirty boson problem in quantum magnets must be taken into account. 

First we note that the singlet - triplet excitation gap A st , proportional to the critical 
field, H c , decreases under high pressure. This was experimentally observed in Ref. [28] for 
the pure spin system TICUCI3. On the other hand it can be argued that the doping acts 
as a chemical pressure, which decreases H c . In fact, since the ionic radius of K + is smaller 
than that of Tl + , a partial substitution of Tl + ions with K + ions produces not only the 
exchange randomness, but also a compression of the crystal lattice. Thus the increase 
of the doping parameter, x, leads to decrease in H c which has indeed been observed 
experimentally [HI |29j [30] . Second, the disorder may increase the triplon effective mass 
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Figure 2. The experimental low temperature magnetization in units of Bohr 
magneton per Cu ion of Tli-zK-rCuCLj for various x in H = 7 T magnetic field 
obtained in Ref . [T4] . 

thereby decreasing the critical temperature T c even when the gap decreases (similar 
effects were observed for helium in porous media [311 [32]). Note that this effect manifests 
itself in different ways. For example, for the mixed compound IPACu(Cl x Bri_ x )3 the 
critical field, H c remains almost unchanged with varying x and then, abruptly becomes 
zero near the Cl-rich phase [33]. In another triplon-BEC compound, Ni(Cli__Br_)2- 
4SC(NH 2 )2, it decreases by a factor of two when x changes from zero to 0.08 [30] 
although the physics of this decrease can be different from that in Tb^Kj-CuCls due 
to the fact that Br atomic radius is larger than the atomic radius of CI. These effects of 
renormalization of the triplon spectrum by disorder can be considered similarly to the 
virtual crystal approximation in the simulations of disorder in solids, where the disorder 
is assumed to lead to a uniform change in the system parameters. The effects of disorder 
such as the appearance of the glassy phase with the density pc and related phenomena 
manifest themselves in addition to these uniform changes. 

The phase diagram of Tb_3.K3.Cud3 in the (H, T) plane was experimentally 
determined in Refs.[__J [29] for various doping x, and the critical field H c was also 
estimated by extrapolation to zero temperature. In the present work the T C (H) 
dependence is given by Eq. (j36l) . We made an attempt to least - square fit our 
parameters m and v by using Eq. ([36]) to describe the experimental phase diagram. 
For simplicity we assume that interparticle interaction is not changed by doping, i.e. 
U = U(R = 0) = 313 K. The parameters obtained by this optimization are presented 
in Table 1. 

Having fixed the input parameters for certain values of x, we are now in the position 
of recalculating the densities as well the magnetization to compare them with the 



The effects of disorder in dimerized quantum magnets in mean field approximations 12 



Table 1. Optimized values of the input parameters of the model: the critical field, H c 
(taken from Ref. |14|). the disorder parameter z/, and the effective mass m for various 
doping x. The critical density, p c , the healing length, A = lA/2m/x, the interparticlc 

1 /3 

distance, d = l/p c and the localization length, L\ oc = d/v arc estimated at H = 7 
T. It is assumed that the doping effects docs not modify the Landc factor g and U . 
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Figure 3. The condensed (a), glassy (b), and supcrfiuid fractions (c) as functions of 
temperature in the YG approximation for various x marked near the plots with the 
input parameters from Tabic 1. The total density of triplons is shown in Fig. 3(d). 



experiment. Figure |3] shows that the doping decreases po an d p s , and increases pc 
as it is expected due to the introduced disorder. Due to change of H c with x, the total 
density of triplons and hence the magnetization, now increases with increasing x in 
accordance with the experiment. One may conclude that the YG approach may well 
describe the effect of disorder to the magnetization, with the additional assumption of 
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an x dependence of the effective mass and the critical field. 

One of the main characteristics of Bose condensed systems is the speed of the 
Bogoliubov mode c, defined here by the Eq. ( lT6j) . which characterizes the propagation of 
collective excitations in the condensate. It is interesting to mention that the magnitude 
of c is large, being only an order of magnitude less than the speed of sound in the 
crystal. This is due to very small trip Ion mass in TICUCI3. Clearly, disorder modifies 
the small-momentum excitation spectrum of the BEC. Estimates of such modification, 
Ac = c — Co, where Co is the speed of the Bogoliubov mode for the pure system, that exist 
in the literature are controversial. For example, perturbative [1] and hydrodynamic [3] 
approaches give Ac > 0, while Ac < was predicted in Refs. jU Ej. In Fig. @] we 
present the corresponding speed for various doping parameters. It can be seen from 
comparison of Fig]4^a) and Fig]4fb) that both MFA approximations considered here 
show the decrease in c with increasing the disorder strength due to the localization 
effects. However, the effect of disorder is small leading to a less than 10 percent decrease 
in the speed of the Bogoliubov mode. 

However, when the spectrum modification by disorder is also taken into account 
by a renormalization of the triplon mass and the gap, as close to the real situation, the 
dispersion of the sound-like mode in fixed magnetic field slightly increases with increasing 
disorder, reaches a maximum and then starts to decrease, (see Figj4b). This behavior is 
caused by interplay between renormalization of the system parameters and localization 
effects. The former tends to increase c, e.g. by increasing /i, and therefore increasing the 
density, while the latter tends to decrease c, e.g. by decreasing the condensed fraction. 
Note that an increase in c with increasing the density was experimentally observed by 
Andrews et al. [34] for the BEC of sodium atoms. This interplay is illustrated in FigJU^c) 
for p and p s . It can be seen that uniform spectrum renormalization first leads to an 
antidepletion effect, increasing these quantities, while the localization effects impair the 
condensation and superfluidity. 

We now consider the question about the existence of a pure Bose glass phase at 
T = 0, which, strictly speaking, should fulfill the following criteria pUl EH EE]: (i) 
gapless in the excitation spectrum, (ii) insulating behavior, i.e. the superfluid fraction, 
p s = 0, (iii) finite compressibility, and (iv) finite density of states. 

In 1970, Tachiki and Yamada [36] have shown that the Heisenberg-like Hamiltonian 
of s = 1/2 dimers can be rewritten as an effective bosonic Hamiltonian. Recently, 
Roscilde and Haas [37] generalized this bosonization procedure taking into account 
disorder and derived a Bose - Hubbard like Hamiltonian usually applied to study "dirty 
bosons" in optical lattices. Applying Fishers ideas [24] we may expect the formation 
of a pure Bose glass phase for doped magnets such as Tb^Kj-CuChj. Although Monte 
Carlo calculations [22| 137] confirmed its existence the experimental confirmation is still 
a matter of debate [201 EH] - We underline here that these Bose glass phases are localized 
out of the BEC phase, i.e. for H < H c . However, in the present work we have been 
mainly concentrating on the region with H > H c where the gapless phase can be realized 
only within the BEC phase. For this case the definition of BG phase may be simplified 
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Figure 4. (a) The speed of sound-like condensate mode at T = 2 K and H = 6 
T as a function of doping parameter x, taking into account solely rcnormalization of 
the system parameters in Table 1 in the YG (solid line) and the HM (dashed line) 
approaches, (b) The same as in Fig. 4(a), now with effects of disorder taken into 
account, (c) The superfluid and condensed fractions (as marked near the plots) in the 
YG (solid lines) and the HM (dashed lines) approximations with the renormalizcd bare 
spectrum parameters. 

as a phase with p 7^ and p s = since the spectrum of the BEC is gapless by itself. 
In searching for such a phase we studied p s and po at T = for various H > H c and x 
and found no pure BG phase with p s = as illustrated in Fig. [5j 

Note also that, as it is seen from Table 1, for moderate values of x considered here, 
the localization length, i.e the mean free path [39] is larger than interparticle distance, 
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Figure 5. The supcrfluid, p s /p (solid lines) and condensed, po/p (dashed lines), 
fractions as a function of the doping parameter x at T — 0, H = 7 T. Upper 
panel corresponds to the YG approximation and lower panel corresponds to the HM 
approximation presented here for comparison. Graphs in plots (a), (c) were calculated 
without bare spectrum rcnormalization, while graphs in plots (b), (d) take into account 
spectrum rcnormalization presented in Table 1. 



L\ oc > d. The physics of the possible BG phase for H < H c fields will be the subject of 
a separate study. 



6. Conclusions 



We reformulated and applied two existing mean field approximations for the "dirty 
boson" problem to study properties of Tb^KzCuCbj quantum magnets. We showed 
that these approaches can qualitatively explain the magnetization data if a certain 
modification of the model parameters similar to the virtual crystal approximation 
usually applied to electron spectrum of alloys is taken into account. 

In fact, random bond effects in mixed dimerized magnetic compounds manifest 
themselves in a dual way: (i) by modification of internal parameters and (ii) by 
localization on random scatterers. Each of these effects could be studied separately 
in an appropriate theory, but they should be taken into account simultaneously for an 
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adequate description of the measured magnetization data. Although the system becomes 
considerably disordered, the Bose-Einstein condensation does not support formation of 
the pure Bose glass phase. The random bonds lead to a nontrivial behavior of the 
sound-like mode speed: when H is fixed and x is experimentally varied, it increases for 
small x, reaches a maximum value and then decreases. While the speed of this mode 
was measured in dilute BEC of sodium atoms a long time ago [31] , it has never been an 
intense focus of research in dimerized quantum magnets [40] . It could be systematically 
studied, for example, by measuring the dispersion relation of the Bogoliubov mode with 
inelastic neutron scattering techniques. 
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Abstract. We study theoretically the effects of disorder on Bose-Einstein 
condensates (BEC) of bosonic triplon quasiparticles in doped dimerized quantum 
magnets. The condensation occurs in a strong enough magnetic field, where the 
concentration of bosons in the random potential is sufficient to form the condensate. 
The effect of doping is partly modeled by a 5 - correlated distribution of impurities, 
which (i) leads to a uniform renormalization of the system parameters and (ii) produces 
disorder in the system with renormalized parameters. This approach can explain 
qualitatively the available magnetization data on the Tli-zK^CuCbj compound taken 
as an example. In addition to the magnetization, we found that the speed of the 
Bogoliubov mode has a maximum as a function of x. No evidence of the pure Bosc 
glass phase has been found in the BEC regime. 



PACS numbers: 75.45. +j, 03.75.Hh, 03.75.Kk 
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1. Introduction 

The effects of disorder on the properties of Bose-Einstein condensates present interesting 
problems both for theoretical and experimental physics [H |2l [31 HI [6]. Disorder 
is important in various systems of real particles such as superfluid 4 He, cold atoms in 
optical lattices, and quasiparticles such as polaritons [7] and excitons [8] • These systems 
are well-suited for experimental studies, however the theory of disordered ensembles of 
interacting bosons is complex and there are essentially no exact solutions even in one 
dimension [9]. To approach this problem, Yukalov and Graham (YG) developed a 
self-consisted stochastic mean field approximation (MFA) [10] for Bose systems with 
arbitrarily strong interparticle repulsion and arbitrary strength of disorder potential. It 
was shown that, in general, the Bose system consists of following coexisting components: 
the condensate fraction, po, the normal fraction p^, the glassy fraction pa, and, in 
addition, can be characterized by the superfluid density p s . In the limit of asymptotically 
weak interactions and disorder the known results, obtained in pioneering work by Huang 
and Meng [11] (HM) are reproduced by the YG theory. An interesting question here 
concerns the problem about the existence of a pure Bose glass (BG) phase, i.e. the 
phase where the condensate fraction is nonzero, while the superfluid fraction is not 
yet present. Note that Ref. [6] introduced an alternative definition of the gapless BG 
phase, having localized short-lived excitations and vanishing superfluid density with 
a continuous transition to the normal phase at finite temperature. Even without 
disorder, the condensate is depleted by particle-particle interactions and temperature. 
The inclusion of random fields depletes the condensate further and, possibly, creates the 
glassy fraction. 

As it was understood recently, a new class of BEC can be provided by spin- 
related quasiparticles in magnetic solids such as intensively pumped magnons [12] or 
triplons in the dimerized quantum magnets in the equilibrium [13] . In the magnets, the 
effect of disorder, which can be produced by admixing other chemical elements, can be 
rather strong to be seen in the physical properties such as the temperature-dependent 
magnetization. The so far most investigated compound showing BEC of triplons is 
TlCuCl 3 . To study the effect of disorder, solid solutions of quantum antiferromagnets 
TICUCI3 and KCuChj, i.e. Tli-^K^CuCls have been experimentally investigated recently 
[HI [151 [16] at low temperatures T. The zero-field ground states of TlCuCl 3 and KCuCl 3 
are spin singlets with excitation gaps A st = 7.1 K and A st = 31.2 K, respectively and 
the magnetic excitations are spin triplets. Triplons arise in magnetic fields H > H c , 
where H c is defined by condition of closing the gap by the Zeeman splitting, that is 
A st = OP'bHc, where g is the electron Lande factor and p,B is the Bohr magneton. In 
the mixture Tli_-K-CuCl3 the induced magnetization M exhibits a cusplike minimum 
at a critical temperature T C (H) for fixed magnetic field H > H c similarly to the parent 
compound, and can be successfully explained in terms of triplon BEC [T7J [18j [19] . 

For a theoretical description it is natural to assume that for weak doping a; < 1 
in the mixed system Tb_a.Ka.Cud3 a small admixture of potassium forms a disorder 
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potential. Consequently, the recently developed theories of "dirty bosons" [T| l_j [TU| [TTj 
can be applied to study the BEC of triplons in Tb^KzCuCls. Here the following 
natural questions arise. For example, what is the correspondence between admixing 
parameter x and the properties of the disorder potential ? What are the experimental 
consequences of the disorder ? Yamada et al. [15] analyzed the electron spin resonance 
spectrum in Tb.ajK^CuCFj and concluded that there is a Bose glass - BEC transition 
near a critical magnetic field. Although this interpretation might need a further analysis 
(see discussion in Refs. [20]) it would be interesting to study the influence of the glassy 
phase, or more exactly, of the glassy fraction pa on the magnetization. Note that even 
the existence of a pure Bose glass phase still is a matter of debate even in theoretical 
approaches. For example, it may be predicted by the approach used by Huang and 
Meng [H] if one extends their formulas from weak disorder to a strong one. On the 
other hand, no pure Bose glass was found in Monte - Carlo simulations [21] for atomic 
gases, but predicted for triplons at T = by Nohadani et al. [22]. 

Here we develop a theory of the disorder effects on the BEC of triplons taking 
Tb_a.Ka.CuO3 as a prototype for studies of specific properties. For example, in atomic 
gases considered in Refs. [T_] [H] , the chemical potential fi is determined self-consistently 
with fixed number of atoms, while in the triplon gas the chemical potential is a given 
external parameter controlled by the applied magnetic field and the number of triplons 
is conserved in the thermodynamic limit. To clarify the terms, we underline that the 
number of magnons may vary but that of triplons may be tuned and kept fixed, which 
makes possible the BEC of the latter. 

The paper is organized as follows. In Sections II and III we outline the YG and 
HM approaches valid only for T < T c and extend it for the triplon system. The shift 
of T c due to disorder and the normal phase properties will be discussed in Section IV. 
Our numerical results will be presented in Section V. Conclusions will summarize the 
results of this work. 



2. Yukalov- Graham approximation for disordered triplons 

In the following we reformulate the Yukalov- Graham approximation to the triplon 
system with arbitrary disorder. The Hamiltonian operator of triplons with contact 
interaction and implemented disorder potential V(r) is given by 

U , , + ^ 



H = d 6 r 



^(r) (k - im + V(r)) ^(r) + - {^(r^r))' 



(1) 



where ip(r) is the bosonic field operator, U is the interparticle interaction strength, and 
K is the kinetic energy operator which defines the bare triplon dispersion e^. Since the 
triplon BEC occurs in solids, we perform integration over the unit cell of the crystal with 
the corresponding momenta defined in the first Brillouin zone. Below we assume that the 
bare spectrum remains coherent in the presence of disorder and consider it as a simple 
isotropic one: = k 2 /2m, where m is the triplon effective mass. The distribution of 
random fields is assumed to be zero - centered, (^(r)) = 0, and the correlation function 
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R(r — r') = (V(r)V(r')). Here and below we adopt the units k B = 1, H = 1, and V = 1 
if not stated otherwise for the unit cell volume. 

To describe Bose condensed system where the global gauge symmetry is broken, 
one employs the Bogoliubov shift: 



V(r) = Vi^W+^iW (2) 

where the condensate density Po( r ) is constant for the homogeneous system, po( r ) = Po- 
Since by the definition the average of ip^(r)ijj(r) is the total number of particles: 



N = d 3 r(^(r)ip(r)) (3) 
Jv 

with the density of triplons per unit cell p = N/V, from the normalization condition 

P = Po + Pi (4) 
one immediately obtains 

\ [ cZ 3 r<v4(r)Vi(r)). (5) 



Therefore the field operator ip\(r) determines the density of uncondensed particles. 

The YG approximation is formulated in representative ensemble formalism, which 
includes two Lagrange multipliers, /i and pi, defined as: 

dpuo dpi 

where Q is the grand thermodynamic potential. It was shown that disorder would not 
change the explicit expressions for chemical potentials, obtained earlier [23] in Hartree- 
Fock-Bogoliubov (HFB) approximation without disorder, 

Po = U{p + Pi Pi = U{p + P\ ~ or), (7) 

1 r « 



where a = — d r (^i(r)tfji(r)) is the anomalous density. The total system chemical 
V Jv 

potential mu, related to the total number of particles as N = —dQ/dfi, is determined 
by 

fip = //ipi + /i po- (8) 

Clearly, when the gauge symmetry is not broken, i.e. po = 0, a = 0, p\ = p, both p 
and pi coincide giving p = pi = 2Up. 

In contrast to homogeneous atomic gases considered in Refs. flDlITT] . where p is fixed 
and p(p) should be calculated as an output parameter, in the trip Ion gas the chemical 
potential is fixed by the external magnetic field, while the density p = p(p) should be 
calculated self consistently. In fact, in a system of triplons p characterizes an additional 
direct contribution to the triplon energy due to the external field H and can be written 
as 

p = gp B H - Am, (9) 
which can be interpreted as a chemical potential of the S z — — 1 triplons. 
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The magnetization is proportional to the triplon density 

M = gn BP (10) 
with p is defined from (jSJ) as 

p = - (/iipi + (jiopo) (11) 
H 

where /i and p,\ are given in © and the densities p , p\ must be calculated self 
consistently. 

It is well known [23] that the disorder field leads to formation of a glassy fraction 
with the density pa- In this approximation each of pi and a are presented as 

Pi = Pn + Pg] (T = o- N + p G (12) 

where p^ and are the normal and anomalous densities without disorder. In the YG 
method, based on HFB approximation, the following explicit relations can be obtained 
HH: 



(Am) 3 / 2 r d s k £k + A 

37T 2 J (27f) £k 



P^v = _ ? + / 7TT3/s(^fcJ— , (13) 



(Am) 3 / 2 [ d 3 k „ . c .l 



with the Bose distribution of Bogoliubov excitations /b(^) = l/(e £fc / T — 1) having the 
dispersion 



^ = vW^fc + 2A. (15) 

For small momentum /c the dispersion is linear, = c/c, and the speed of the Bogoliubov 
mode 

\/A , . 

Vm 

The self energy A is determined formally by the same equation as in the case when 
the disorder is neglected, 

A = U(po + a) = U(p- Pn + o- n ). (17) 

The contribution from the disorder potential is hidden in the density of the glassy 
fraction 

p G = ^ J d 3 r((M*)Mr))) (18) 

where the double angle brackets mean the stochastic average. In general, the calculation 
of pa is rather complicated, however, for the 5 - correlated disorder i.e. for the white 
noise, 

((V(r)V(r'))) = R5(r-r'), (19) 
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equation (ITS!) is simplified as [TU] 

_ Ro(p ~ Pn) , , 

PG ~ R + 7(l-R )^- (20) 

The density of condensed fraction can be found by inserting (12) and (20) into the 
normalization condition fl4]). The result is 

7(1 -Rof/\ P - PN ) 
p0 ~ R + 7(l-R )^ ■ (21) 

In Eqs. (1201) and (1211) we introduced the dimensionless parameter Rq as 

7Rm 2 , . 

Ro = — r =. (22) 
47TV mA 

One can see from Eqs. ( !20|) and ( T2~TT) that the glassy fraction is proportional to the 
condensed one, 

pG ~ 7{i-R yn- (23) 

The system of Eqs. (JTJ) , ([8]) , (JT3l)- (fT9l) are the basic of YG approximation. 

An interesting quantity, crucial for determining the Bose glass phase, is the 
superfluid density, p s . In general it is defined as a partial density appearing as a response 
to a velocity boost 

ft = S57S5^ A > (24) 

where P v is the total momentum of the system, dependent on the macroscopic velocity 
v. Referring the reader to original papers JTOj [11] we bring below analytical expression 
obtained there for p s in the case of white noise random potential 

Ps =,- — -—, (25) 

i r k 2 d 3 k 

lm~J (27r) 3 sinh 2 (^ fc /2T) 
Note that YG approach is valid for arbitrary strength of the interaction potential 
U, and for arbitrary strong disorder. For the weak interactions it leads to pioneering 
Huang- Meng approach which will be extended to the "dirty triplons" in the next 
section. 



Qn = — I 77^—727^17^- ( 26 ) 



3. Huang-Meng approximation 

For completeness, we present here the results for the Huang-Meng approach, based on 
the so called Hartree Fock Popov (HFP) approximation which has been widely applied 
in the literature to describe the BEC of triplons [171 [19]. The basic equations of this 
approach can be obtained by neglecting the anomalous density a, which leads naturally 
to the single chemical potential p, = fio = P>i- Namely, one finds from (JTJ), (JH]) and ( TT71) 



A = Upo, p = U(p + Pl ). 



(27) 
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From these equations and (1121) one obtains following main equations for the self energy 
A: 

A = p-2U(p N + p G ), (28) 

where p^ is formally given in (I13p . and po is determined by the first equation in fl27|) . 
The glassy fraction can be obtained from (1201) in the linear approximation by R assuming 
weakness of interparticle interaction [10] E] 

where a s = Um/An is the s - wave scattering length. Inserting (|29|) into (|28l) we can 
rewrite the former as 

m 2 R\ r A . . 

A = H-2Up N - n L_ . 30 
27r-v/m 



To evaluate the densities one has to solve nonlinear algebraic equation ( 130]) . where p^ 
is given formally by (!I3|) . with respect to A. Next, by inserting the result into ( |27l) 
and ( 129]) one obtains the density of condensed triplons p an d the glassy fraction p G , 
respectively. The total density can be evaluated then by the normalization condition 
P = Po + Pn + Pg- Equations (1231) . ( 126]) for the superfluid density are formally the same 
in both approximations. 

4. The shift of the critical temperature due to disorder and the T > T c 
regime 

It is well known that the critical temperature of BEC, T c for an ideal gas is given by: 

where p c is the total density of triplons near the critical temperature of BEC for pure 

system, 

p c = n/2U, (32) 

with C{x) being the Riemann function. Equation (132]) directly follows, from Eqs. (J7J), 
(JHJ) or (128]) by setting p N = p and p = p G = 0. 

Clearly, any type of interaction is expected to modify T c . In general, these 
modifications are related to the interparticle interactions as well as to the disorder 
potential. Both approaches, considered here give a zero shift due to the boson-boson 
repulsion. However the shift due to the S— correlated disorder (fT9]) . AT C = T c — T c ° is 
given as [TQl EH] 

AT r 2u 



T c ° 9tt' 
where the dimensionless disorder parameter v 

1 



1/3 T ' 
Pc Moc 



(33) 



(34) 
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c 



is introduced with the localization length 

Lloc = WR- (35) 
For practical calculations we rewrite T c in Eq. (|33|) . which is in a good agreement with 
perturbative estimates [2] as well as with Monte Carlo simulations [26], as an explicit 
function of effective mass m, the interaction strength U, critical magnetic field H c , 
disorder parameter u, and external field H as follows: 

9tt - 2u ( V2g/2 B (H - H ( " 

Now we proceed to consider the triplon density in the normal state in the T — T c ^> AT, 
temperature range. The dirty bosons in the normal phase where the gauge symmetry 
is not broken, are yet poorly studied. For R = with p = pc = o = the triplon gas 
behaves like an "ideal gas" with an effective chemical potential p e ff , and the density 

d 3 k 1 
(27r) 3 exp((e fc -/i cfr )/T) - 1 
Although p, e ff is not accurately known it depends in general, on p, as well as on R. For 
the pure case MFA [27] gives p e s(R = 0) = p — 2Up. The contribution from the disorder 
potential has been studied neither in YM nor in HM approaches. Therefore, to make 
the calculations self consistently, we have to use 

f d 3 k 1 
P(T >Tc) -J (27f exp ((e k - p - 2U P )/T) - 1 ' (38) 
which yields the density p as a solution of the nonlinear equation (138j) . 



p(T >T C )= I „ Wfm r - (37) 



5. Results and discussions 

In the calculations below, the energies are measured in Kelvin, the mass in K _1 , the 
densities are dimensionless and the Bohr magneton is fiB = 0.671668 K/T. As to the 
strength of disorder potential R, it has units K -2 while the disorder parameter u, defined 
in Eq. (1511) is a number supposed to be less than one, v < 1. As a material parameter, 
we use mean dimer-dimer distance in TlCuChj r^d = 0.79 nm [19]. 

To perform numerical calculations in the YG approximation, assuming that fi, U, 
m, and R are given parameters, we use following strategy, (i) By inserting (J7J, (fl2]) . 
(120]) and (T2T1) into (ITT]) we obtain quadratic algebraic equation with respect to p and 
solve it analytically, (ii) By using this p(fi,R,A) and ( IT3]) . ( IT4l) in ( IT71) we solve the 
latter numerically with respect to A, and (iii) by inserting this A back into p(p, R, A) 
we find the magnetization from ( ITOl and evaluate other densities like po and Pg from 
(|2Q) and d2U>. 

In Figure UJ we present as an example the total triplon density p(T) for a clean and 
strongly disordered [y = 0.45, see Eq.( l34l) ) TlCuCFj, obtained in the YG approximation 
assuming that the total effect of the doping leads only to the randomness in the triplon 
subsystem. 
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Figure 1. The total triplon density as a function of temperature in the YG 
approximation for two values of v. Here the following set of parameters m=0.0204 
K _1 , A st = 7.3 K,U = 313 K, and g = 2.06 [Hj valid for TlCuCl 3 is used. 

The calculation of other quantities using the same assumption shows that the 
disorder leads to a decrease in the condensed and superfluid fractions, thereby increasing 
the glassy one. This tendency is quite natural, since the localization effects prevent 
particles from going into BEC. However, the increase in pc is so weak that along 
with p the total number of triplons p is also decreased with increasing the strength 
of disorder potential R. Bearing in mind that p is proportional to the magnetization M, 
and v is assumed to be approximately proportional to x, and comparing FigJT]with the 
experimental data illustrated in Fig J2] one may conclude that the agreement between the 
theory and the experiment is unsatisfactory since the main features of the experimental 
results are not reproduced there. As it is seen in Figj2]the disorder leads to an increase 
in the magnetization and, hence, in the total triplon density. This is accompanied 
by the decrease in the transition temperature. We therefore conclude that while the 
triplon gas can be considered similarly to atomic gases for which the considered mean- 
field approximations were developed, some further additional specific material - related 
properties of the dirty boson problem in quantum magnets must be taken into account. 

First we note that the singlet - triplet excitation gap A st , proportional to the critical 
field, H c , decreases under high pressure. This was experimentally observed in Ref. [28] for 
the pure spin system TICUCI3. On the other hand it can be argued that the doping acts 
as a chemical pressure, which decreases H c . In fact, since the ionic radius of K + is smaller 
than that of Tl + , a partial substitution of Tl + ions with K + ions produces not only the 
exchange randomness, but also a compression of the crystal lattice. Thus the increase 
of the doping parameter, x, leads to decrease in H c which has indeed been observed 
experimentally [HI |29j [30] . Second, the disorder may increase the triplon effective mass 
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Figure 2. The experimental low temperature magnetization in units of Bohr 
magneton per Cu ion of Tli-zK-rCuCLj for various x in H = 7 T magnetic field 
obtained in Ref . [T4] . 

thereby decreasing the critical temperature T c even when the gap decreases (similar 
effects were observed for helium in porous media [311 [32]). Note that this effect manifests 
itself in different ways. For example, for the mixed compound IPACu(Cl x Bri_ x )3 the 
critical field, H c remains almost unchanged with varying x and then, abruptly becomes 
zero near the Cl-rich phase [33]. In another triplon-BEC compound, Ni(Cli__Br_)2- 
4SC(NH 2 )2, it decreases by a factor of two when x changes from zero to 0.08 [30] 
although the physics of this decrease can be different from that in Tb^Kj-CuCls due 
to the fact that Br atomic radius is larger than the atomic radius of CI. These effects of 
renormalization of the triplon spectrum by disorder can be considered similarly to the 
virtual crystal approximation in the simulations of disorder in solids, where the disorder 
is assumed to lead to a uniform change in the system parameters. The effects of disorder 
such as the appearance of the glassy phase with the density pc and related phenomena 
manifest themselves in addition to these uniform changes. 

The phase diagram of Tb_3.K3.Cud3 in the (H, T) plane was experimentally 
determined in Refs.[__J [29] for various doping x, and the critical field H c was also 
estimated by extrapolation to zero temperature. In the present work the T C (H) 
dependence is given by Eq. (j36l) . We made an attempt to least - square fit our 
parameters m and v by using Eq. ([36]) to describe the experimental phase diagram. 
For simplicity we assume that interparticle interaction is not changed by doping, i.e. 
U = U(R = 0) = 313 K. The parameters obtained by this optimization are presented 
in Table 1. 

Having fixed the input parameters for certain values of x, we are now in the position 
of recalculating the densities as well the magnetization to compare them with the 
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Table 1. Optimized values of the input parameters of the model: the critical field, H c 
(taken from Ref. |14|). the disorder parameter z/, and the effective mass m for various 
doping x. The critical density, p c , the healing length, A = lA/2m/x, the interparticlc 

1 /3 

distance, d = l/p c and the localization length, L\ oc = d/v arc estimated at H = 7 
T. It is assumed that the doping effects docs not modify the Landc factor g and U . 
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Figure 3. The condensed (a), glassy (b), and supcrfiuid fractions (c) as functions of 
temperature in the YG approximation for various x marked near the plots with the 
input parameters from Tabic 1. The total density of triplons is shown in Fig. 3(d). 



experiment. Figure |3] shows that the doping decreases po an d p s , and increases pc 
as it is expected due to the introduced disorder. Due to change of H c with x, the total 
density of triplons and hence the magnetization, now increases with increasing x in 
accordance with the experiment. One may conclude that the YG approach may well 
describe the effect of disorder to the magnetization, with the additional assumption of 



The effects of disorder in dimerized quantum magnets in mean field approximations 13 

an x dependence of the effective mass and the critical field. 

One of the main characteristics of Bose condensed systems is the speed of the 
Bogoliubov mode c, defined here by the Eq. ( lT6j) . which characterizes the propagation of 
collective excitations in the condensate. It is interesting to mention that the magnitude 
of c is large, being only an order of magnitude less than the speed of sound in the 
crystal. This is due to very small trip Ion mass in TICUCI3. Clearly, disorder modifies 
the small-momentum excitation spectrum of the BEC. Estimates of such modification, 
Ac = c — Co, where Co is the speed of the Bogoliubov mode for the pure system, that exist 
in the literature are controversial. For example, perturbative [1] and hydrodynamic [3] 
approaches give Ac > 0, while Ac < was predicted in Refs. jU Ej. In Fig. @] we 
present the corresponding speed for various doping parameters. It can be seen from 
comparison of Fig]4^a) and Fig]4fb) that both MFA approximations considered here 
show the decrease in c with increasing the disorder strength due to the localization 
effects. However, the effect of disorder is small leading to a less than 10 percent decrease 
in the speed of the Bogoliubov mode. 

However, when the spectrum modification by disorder is also taken into account 
by a renormalization of the triplon mass and the gap, as close to the real situation, the 
dispersion of the sound-like mode in fixed magnetic field slightly increases with increasing 
disorder, reaches a maximum and then starts to decrease, (see Figj4b). This behavior is 
caused by interplay between renormalization of the system parameters and localization 
effects. The former tends to increase c, e.g. by increasing /i, and therefore increasing the 
density, while the latter tends to decrease c, e.g. by decreasing the condensed fraction. 
Note that an increase in c with increasing the density was experimentally observed by 
Andrews et al. [34] for the BEC of sodium atoms. This interplay is illustrated in FigJU^c) 
for p and p s . It can be seen that uniform spectrum renormalization first leads to an 
antidepletion effect, increasing these quantities, while the localization effects impair the 
condensation and superfluidity. 

We now consider the question about the existence of a pure Bose glass phase at 
T = 0, which, strictly speaking, should fulfill the following criteria pUl EH EE]: (i) 
gapless in the excitation spectrum, (ii) insulating behavior, i.e. the superfluid fraction, 
p s = 0, (iii) finite compressibility, and (iv) finite density of states. 

In 1970, Tachiki and Yamada [36] have shown that the Heisenberg-like Hamiltonian 
of s = 1/2 dimers can be rewritten as an effective bosonic Hamiltonian. Recently, 
Roscilde and Haas [37] generalized this bosonization procedure taking into account 
disorder and derived a Bose - Hubbard like Hamiltonian usually applied to study "dirty 
bosons" in optical lattices. Applying Fishers ideas [24] we may expect the formation 
of a pure Bose glass phase for doped magnets such as Tb^Kj-CuChj. Although Monte 
Carlo calculations [22| 137] confirmed its existence the experimental confirmation is still 
a matter of debate [201 EH] - We underline here that these Bose glass phases are localized 
out of the BEC phase, i.e. for H < H c . However, in the present work we have been 
mainly concentrating on the region with H > H c where the gapless phase can be realized 
only within the BEC phase. For this case the definition of BG phase may be simplified 
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Figure 4. (a) The speed of sound-like condensate mode at T = 2 K and H = 6 
T as a function of doping parameter x, taking into account solely rcnormalization of 
the system parameters in Table 1 in the YG (solid line) and the HM (dashed line) 
approaches, (b) The same as in Fig. 4(a), now with effects of disorder taken into 
account, (c) The superfluid and condensed fractions (as marked near the plots) in the 
YG (solid lines) and the HM (dashed lines) approximations with the renormalizcd bare 
spectrum parameters. 

as a phase with p 7^ and p s = since the spectrum of the BEC is gapless by itself. 
In searching for such a phase we studied p s and po at T = for various H > H c and x 
and found no pure BG phase with p s = as illustrated in Fig. [5j 

Note also that, as it is seen from Table 1, for moderate values of x considered here, 
the localization length, i.e the mean free path [39] is larger than interparticle distance, 
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Figure 5. The supcrfluid, p s /p (solid lines) and condensed, po/p (dashed lines), 
fractions as a function of the doping parameter x at T — 0, H = 7 T. Upper 
panel corresponds to the YG approximation and lower panel corresponds to the HM 
approximation presented here for comparison. Graphs in plots (a), (c) were calculated 
without bare spectrum rcnormalization, while graphs in plots (b), (d) take into account 
spectrum rcnormalization presented in Table 1. 



L\ oc > d. The physics of the possible BG phase for H < H c fields will be the subject of 
a separate study. 



6. Conclusions 



We reformulated and applied two existing mean field approximations for the "dirty 
boson" problem to study properties of Tb^KzCuCbj quantum magnets. We showed 
that these approaches can qualitatively explain the magnetization data if a certain 
modification of the model parameters similar to the virtual crystal approximation 
usually applied to electron spectrum of alloys is taken into account. 

In fact, random bond effects in mixed dimerized magnetic compounds manifest 
themselves in a dual way: (i) by modification of internal parameters and (ii) by 
localization on random scatterers. Each of these effects could be studied separately 
in an appropriate theory, but they should be taken into account simultaneously for an 
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adequate description of the measured magnetization data. Although the system becomes 
considerably disordered, the Bose-Einstein condensation does not support formation of 
the pure Bose glass phase. The random bonds lead to a nontrivial behavior of the 
sound-like mode speed: when H is fixed and x is experimentally varied, it increases for 
small x, reaches a maximum value and then decreases. While the speed of this mode 
was measured in dilute BEC of sodium atoms a long time ago [31] , it has never been an 
intense focus of research in dimerized quantum magnets [40] . It could be systematically 
studied, for example, by measuring the dispersion relation of the Bogoliubov mode with 
inelastic neutron scattering techniques. 
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